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Abstract 


We apply heavy-quark effective theory to separate long- and short-distance 
effects of heavy quarks in lattice gauge theory. In this approach, the inverse 
heavy-quark mass and the lattice spacing are treated as short distances, and 
their effects are lumped into short-distance coefficients. We show how to 
use this formalism to match lattice gauge theory to continuum QCD, order 
by order in the heavy-quark expansion. In this paper, we focus on heavy- 
light currents. In particular, we obtain one-loop results for the matching 
factors of lattice currents, needed for heavy-quark phenomenology, such as 
the calculation of heavy-light decay constants, and heavy-to-light transition 
form factors. Results for the Brodsky-Lepage-Mackenzie scale q* are also 
given. 
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I. INTRODUCTION 


A key ingredient in flavor physics is the calculation of hadronic matrix elements of the 
electroweak Hamiltonian. For example, one would like to calculate, from first principles, 
quantities such as leptonic decay constants, semi-leptonic form factors, and the amplitudes 
for neutral-meson mixing. Numerical calculations with lattice QCD offer a way to obtain 
these quantities, eventually with well-controlled estimates of the numerical uncertainties [1]. 

The properties of B and D mesons are especially interesting, but the relatively large b 
and c quark masses make it difficult, with today’s computers, to carry out lattice calculations 
in the limit mga —0 for which lattice QCD was first developed. (Here rriQ is the 6 or c 
quark mass, and a is the lattice spacing.) One can, however, use the simplifying features 
of the heavy-quark limit of QCD to make lattice calculations tractable. As mg is increased 
far above the typical scale of the wave function, Aqcd, the hadrons’ wave functions depend 
less and less on mg. As mg oo the wave functions become flavor and spin symmetric [2]. 
For quarkonia similar simplifications occur, including spin symmetry [3]. 

In this paper we construct vector and axial vector currents with one quark heavy and the 
other light. These currents are needed to obtain the decay constants of heavy-light mesons, 
and the form factors for decays of the form H Llui, where H is a charmed or 6-flavored 
hadron {e.g., B, D; A^, Ac), decaying to a light hadron L {e.g., tt, K, p; p, etc.) and a lepton I 
and its neutrino z/;. In particular, we provide a way to treat radiative and power corrections 
consistently. This paper is a sequel to Ref. [4], which focussed on power corrections. Here we 
discuss the case of heavy-light bilinears in detail, and we compute explicitly the matching 
factors for the currents introduced in the “Fermilab” formalism [5]. Heavy-heavy bilinears 
are considered in a companion paper [6]. 

To interpret lattice calculations when mga 1, it is convenient to describe cutoff effects 
with the Symanzik local effective Lagrangian (LE£) and expand the LE£’s short-distance 
coefficients in powers of mga [7-10]. When mga 1, however, one should realize that it is 
not lattice gauge theory that breaks down but rather the Symanzik description, especially its 
expansion in mga. If mga is large because mg 3> Aqcd, then the simplifying features of the 
heavy-quark limit provide an alternative. Instead of matching lattice gauge theory directly to 
continuum QCD, one can match to the heavy-quark effective theory (HQET) or, for quarko¬ 
nia, to non-relativistic QCD (NRQCD). In this approach, the inverse heavy-quark mass and 
the lattice spacing are both treated as short distances, and a simple picture arises, in which 
heavy-quark discretization effects are lumped into short-distance coefficients. Heavy-quark 
cutoff effects are systematically reducible, by adjusting the heavy-quark expansion for lattice 
gauge theory to agree term-by-term with continuum QCD. 

Such application of HQET to lattice QCD was started in Ref. [4], building on Ref. [5]. 
In this paper we extend the formalism to heavy-light currents. We use the heavy-quark 
expansion, as generated by HQET, to derive matching conditions, which are valid for all 
mga and to all orders in the gauge coupling. Our derivation is explicit for dimension-four 
currents, which is the next-to-leading dimension, but generalization to higher-dimension 
operators should be clear. 

We also present explicit results for the one-loop radiative corrections to the normalization 
of the current. These calculations show that the temporal and spatial components of the cur¬ 
rent do not have the same radiative corrections. This feature has been found already [11,12], 
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and the HQET formalism shows why it arises. In deriving these results we have found a 
compact way of arranging the Dirac algebra, which may be useful for calculations with 
actions, such as highly improved actions, that are not considered here. 

As expected, the coefficients have a strong mass dependence. Most of this dependence 
can be handled non-perturbatively [13-15]. For equal mass, it is simple to normalize the 
temporal vector current, for all masses. One can then form ratios of renormalization factors, 
from which the dominant mass dependence drops out. Results for these combinations are 
also presented, in Sec. IV. 

Our one-loop results extend those of Ref. [12], which considered heavy-light currents 
with the Sheikholeslami-Wohlert (SW) action [16] for Wilson fermions [17] and also with 
non-relativistic QCD (NRQCD). Results for the Wilson action [17] have been obtained first 
by Kuramashi [11]. In Refs. [11,12] a term in the currents, the so-called rotation term [18,5], 
which is needed for tree-level improvement at order l/rug, was omitted. Here we include the 
rotation, obtaining the algebraic expression of the Feynman diagrams for the full Fermilab 
action. We present numerical results for the Wilson action (without rotation) and the SW 
action (with and without rotation). These results are appropriate for recent calculations 
of decay constants [19-23], which used the radiative corrections calculated in Refs. [11,12]. 
Our new results have been used in a recent calculation of the form factors for the decays 
B —> ttIvi and D —> ttIvi [15]. We also have obtained results for the Fermilab action on 
anisotropic lattices [24]. 

Our formalism should be useful for computing matching factors (beyond one-loop) also in 
lattice NRQCD [25]. Applied to the static limit [26], it generalizes the formalism of Eichten 
and Hill [27]. At one-loop order, similar methods have been developed to calculate the heavy- 
light matching coefficients for lattice NRQCD [28,29]. As in the Symanzik program [7-9], 
the advantage of introducing a continuum effective field theory is that the formalism pro¬ 
vides a clear dehnition of the matching coefficients at every order in perturbation theory 
(in the gauge coupling). Indeed, it may also provide a foundation for a non-perturbative 
improvement program. 

This paper is organized as follows. Section H discusses three ways to separate long and 
short distance physics with (continuum) effective held theories. The hrst is Symanzik’s 
description of lattice spacing effects; we also discuss its breakdown when mqa 1. The 
second is the HQET description of heavy quarks, applied to continuum QCD. The third is 
the HQET description of heavy quarks on the lattice, which applies when mq 3> Aqcd, for 
all mqa. In particular, we obtain a dehnition of the matching factors for the vector and 
axial-vector heavy-light currents. Section H also shows how the HQET matching procedure 
is related to the Symanzik procedure in the regime where both apply. Then, the Fermilab 
action is reviewed in Sec. HI, and in Sec. IV we present one-loop results for the matching 
factors. Some concluding remarks are made in Sec. V. Three appendices contain details of 
the one-loop calculation, including an outline of a method to obtain compact expressions, 
and explicit results for the one-loop Feynman integrands for the renormalization factors with 
the full Fermilab action. 

Instead of printing tables of the numerical results in Sec. IV, we are making a suite of 
programs freely available [30]. This suite includes programs for the heavy-heavy currents 
treated in our companion paper [6]. 
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II. MATCHING TO CONTINUUM FIELD THEORIES 


In this section we discuss how to interpret the physical content of lattice field theories by 
matching to continuum field theories. First, the standard Symanzik formalism for describing 
cutoff effects is reviewed, and we recall how this description breaks down for heavy quarks. 
After reviewing the HQET description of (continuum) QCD, we adapt HQET to describe 
lattice gauge theory. Comparison of the two then yields a matching procedure that is valid 
whenever uiq 3> Aqcd, and for all rriga. In the limit mqa 1 both the HQET and the 
Symanzik descriptions should apply, so we are able to derive relations between the some of 
the matching coefficients. 


A. Symanzik Formalism 

The customary way to define matching factors for lattice gauge theory is to apply 
Symanzik’s formalism. Then the short-distance lattice artifacts are described by a local 
effective Lagrangian (LE£) and local effective operators. For the Lagrangian of any lattice 
field theory one can write [7,8] 


-^lat — (2.1) 

where the symbol = can be read “has the same on-shell matrix elements as”. The left-hand 
side is a lattice field theory, and the right-hand side is a continuum field theory, whose 
ultraviolet behavior is regulated and renormalized completely separately from the lattice of 
the left-hand side. The LE£ is the Lagrangian of the corresponding continuum field theory, 
plus extra terms to describe discretization effects. For lattice QCD 

C-Syra = >CqcD + ^,1-, (2.2) 

where Dqcd is the renormalized, continuum QCD Lagrangian. We focus on the quarks, so 
for our purposes 


>^^QCD = -q{p + TUq) q. (2.3) 

Lattice artifacts are described by higher-dimension operators. 

Cl = aK^.Fqia^^F^^^'q H-, (2.4) 

where a is the lattice spacing and K^.p is a short-distance coefficient that depends on details 
of the lattice action [9]. The lattice artifacts in Cj can be treated as a perturbation. In this 
way a series can be developed, with matrix elements in the (continuum) eigenstates of Dqcd- 

Equation (2.2) omits dimension-five operators of the form qR{p+mq)q or q{—p+mq)Rq, for 
arbitrary R, which make no contribution to on-shell matrix elements, owing to the equations 
of motion implied by Eq. (2.3). 

The vector and axial vector currents can be described in a similar way. Consider, for 
example, the flavor-changing transition s —> u. Then one may write [9] 
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- aKvd^uaf^'^s + ■■■, (2.5) 

+ aKAd^ui'^^s + ■ ■ •, ( 2 . 6 ) 

where 

= uiYs, (2.7) 

= mY75S, ( 2 . 8 ) 


are the vector and axial vector currents in QCD. Further dimension-four operators are 
omitted, because they are linear combinations of those listed and others that vanish by the 
equations of motion. Like the terms of dimension hve and higher in Cj, the dimension-four 
currents can be treated as perturbations. Matrix elements of and then give 

those of continuum QCD, at least in the limit a —0. 

The short-distance coefficients— K^.p, Kj, and Zj ( J = V, A) —are, in general, functions 
of the gauge coupling and the quark masses (in lattice units), and they depend on the 
renormalization scheme of the LE£. For m^a 1 (g = m, d, s), it is consistent and 
satisfactory to replace K^.p and Kj with their values at rUga = 0, and to replace the Zj 
with the hrst two terms of the Taylor expansion around m^a = 0. For example, with Wilson 
fermions [17,16] and conventional bilinears for the lattice currents, one hnds Ky^ = = 0, 

and 

= i(l - csw) + C>(ma), (2.9) 

zP = Z^A = 1 + + rns)a + O(m^a^), (2-10) 

where the superscript “[0]” denotes the tree level, and csw is the clover coupling of the 
SW action [16] (cf. Sec. III). Moreover, in the hands of the Alpha Collaboration [9,10], 
Eqs. (2.1)-(2.6) are the foundation of a non-perturbative procedure for adjusting K„.p, Ky, 
and Ka to be of order aMp, where Mp is a (light) hadronic mass scale, and also for computing 
Zy and Za non-perturbatively (through order Mpo). Then all lattice artifacts in the mass 
spectrum, decay constants, and form factors are of order a^. 

Eor a heavy flavor Q, however, it is not practical to keep mqa small enough so that this 
program straightforwardly applies. Recent work that uses the fully 0(a)-improved action 
and currents has chosen the heavy-quark mass mqa to be as large as 0.7 or so. Thus, {mqo)^ 
is not small,^ and one should check whether contributions of order {mqaY are under control. 
Indeed, if one keeps the full mass dependence in the coefficients, one hnds that the simple 
description of Eqs. (2.2)-(2.6) breaks down. The relation between energy and momentum 
becomes [5,31] 

E‘^{p) = + + Oip^a^), ( 2 . 11 ) 

1712 

where, for the Wilson and SW actions. 


^Also, the lowest chosen values of mq are around 1 GeV, which may be too small to be considered 
“heavy”. 
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( 2 . 12 ) 

(2.13) 


m^i^a = ln(l + rrioa), 

1 _ 2 1 

rnoa{2 + moo) 1 + moa ’ 

and tuq is the bare lattice mass. Generalizations valid at every order in perturbation theory 
also have been derived [31]. In a similar vein, the spatial and temporal components of 
the currents no longer take the same matching coefficients, as shown by explicit one-loop 
calculations [ 11 , 12 ]. 

The energy-momentum relation in Eq. (2.11) is obtained for pa 1 but rriga ^ 1. It 
can be described by modifying the standard LE£ to 



Aat = -q liDi -F W—7 ■ D + mi\q + 


(2.14) 


that is, temporal and spatial directions must be treated asymmetrically in the dimension- 
four Lagrangian, and also in the higher-dimension terms Cj. From the tree-level formulas, 
Eqs. (2.12) and (2.13), 


nil 

m2 


1 — ‘^m\a^ + + • • •, 


(2.15) 


SO one sees that the deviation from the standard description is of order (ma)^. (At the one- 
loop order [31], and at every order in Eq. (2.15) still has no term linear in ma.) One can 
arrive at Eqs. (2.14) and (2.15) also by starting with Eq. (2.2), including higher-dimension 
terms, and eliminating 74 /I 4 and etc., by applying the equations of motion. 

In any case, deviations of mi/m 2 —and similar ratios—from 1 are present in lattice 
calculations. With the Wilson or SW actions 1 — mi/m 2 , for example, is 10 percent or 
greater for mia > 0.6. Although this numerical estimate is made at the tree level, it is 
implausible that radiative corrections or bound-state effects could wash the error away. In 
summary, the description of Eqs. (2.2)-(2.6) is no longer accurate when mga ^ 1. 

There are several possible remedies. One is to do numerical calculations with a so small 
that, even for the b quark, m^a 1. Despite the exponential growth in computer power, 
this remedy will not be available for many years. Another remedy is to add a parameter 
to the lattice action, which can be tuned to set mi = m 2 [5]. An example of this is an 
action with two hopping parameters. Then, the continuum description can again take the 
form in Eq. (2.2), starting with the continuum Dqcd, although it is still useful to describe 
the higher-dimension terms asymmetrically. A third remedy is to realize that it is the 
description, rather than the underlying lattice gauge theory, that has broken down. Since 
lattice gauge theory with Wilson fermions has a well-behaved heavy-quark limit [5], it is 
possible to use heavy-quark effective theory (HQET) or NRQCD to describe short-distance 
effects, including the lattice artifacts of the heavy quark [4]. This last remedy is explained 
in detail in Sec. IIC, where we show also how all three strategies are connected. 


B. HQET description of QCD 

The breakdown of the standard Symanzik description of cutoff effects for Wilson fermions 
arises because the kinematics of heavy hadron decays single out a vector, namely, the heavy 
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hadron velocity. But, since the heavy-quark mass is also much larger than the spatial 
momenta of the problem, the dynamics simplify. In continuum QCD, this has led to the 
development of the effective held theories HQET [26,27,34-36] and NRQCD [3,25]. These 
two effective theories are useful for generating an expansion in p/mq. They share a common 
effective Lagrangian, but the power in p/mq assigned to any given operator is not necessarily 
the same. In HQET the power can be deduced immediately from the dimension, whereas in 
NRQCD it is deduced by counting powers of the relative velocity of the QQ system. The 
discussion in this paper will follow the counting of HQET, but the logic could be repeated 
with the counting of NRQCD. 

Our aim is to show, for the case of heavy-light currents, how to use HQET to extend the 
standard Symanzik program into the region where mqa is no longer small. This program 
was started in Ref. [4], building on Ref. [5]. The formalism holds for all mqa, but, like the 
usual HQET, it requires 


mQ>p, Aqcd- ( 2 . 16 ) 

First, in this subsection, we recall the HQET description of continuum QCD, paralleling the 
discussion in Sec. HA. Then, in Sec. HC, we explain what changes are needed to describe 
the cutoff effects of lattice NRQCD and of lattice gauge theory with Wilson fermions. 

The HQET conventions are the same as those given Sec. HI of Ref. [4]. The velocity 
needed to construct HQET is v. The fourth Euclidean component ^4 = so in the rest 
frame v = {i, 0 ). The metric is taken to be diag(±l, 1,1,1), with the upper (lower) sign for 
Euclidean (Minkowski) spacetime. In either case, = —1. The heavy quark held is called 
and it satishes the constraint 1(1 — or 

= ihy, = ih^. (2-17) 

Physically Eq. (2.17) means that describes only quarks, but not anti-quarks. The tensor 
projects onto components orthogonal to v. For a vector p, the component 
orthogonal to n is p[[ = = p^ -f- n^(n-p); in the rest frame, these are the spatial 

components. 

HQET describes the dynamics of heavy-light bound states with an effective Lagrangian 
built from h^. So, for these states, one can say 

l^QCD = l^HQET, ( 2 . 18 ) 


where 


-Chqet = + ■ ■ ■ • (2.19) 

For HQET contains terms of dimension 4 -|- s. Note that the ultraviolet regulator 
and renormalization scheme of the two sides of Eq. (2.18) need not be the same, although 
dimensional regularization and the MS scheme are usually used for both. 

For this paper it is enough to consider the hrst two terms, and The leading, 
dimension-four term is 


= hv{iv ■ D — m)hy. 


( 2 . 20 ) 
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The choice of v is somewhat arbitrary. If v is close to the heavy quark’s velocity,^ then C^^'> 
is a good starting point for the heavy-quark expansion, which treats the higher-dimension 
operators as small. The most practical choice is the containing hadron’s velocity. 

The mass term in is often omitted. By heavy-quark symmetry, it has an effect 
neither on bound-state wave functions nor, consequently, on matrix elements. It does affect 
the mass spectrum, but only additively. Including the mass obscures the heavy-quark flavor 
symmetry, but only slightly [4]. When the mass term is included, higher-dimension operators 
are constructed with — imv^ [32]. To describe on-shell matrix elements one 

may omit operators that vanish by the equation of motion, —iv ■ T>h^ = 0, derived from 
Eq. (2.20). In practice, therefore, higher-dimension operators are constructed fromT>[j[ = 
and 

The dimension-hve interactions are 


C^^'> =C202 +Cs0s, 

(2.21) 

where C 2 and Cs are short-distance coefficients, and 


O 2 = 

(2.22) 

Os = hvSaf3B°'^K, 

(2.23) 

with Saf 3 = —ia 0 , 13/2 and . 



In Eq. (2.20) one should think of the quark mass m as a short-distance coefficient. By 
reparametrization invariance [37], the same mass appears in the denominator of the kinetic 
energy €2021 namely. 



If operator insertions in HQET are renormalized with a minimal subtraction in dimensional 
regularization, then m is the (perturbative) pole mass. With other ultraviolet regulators, the 
operator and the mass m could become /^-dependent. Even in mass-independent schemes, 
the chromomagnetic operator Os depends on the renormalization point fi of the HQET, and 
that dependence is canceled by 




2 m 


(2.25) 


with 2 m appearing so that zs is unity at the tree level. 

The description of electroweak flavor-changing operators proceeds along the same lines. 
The flavor-changing vector current for a b ^ q transition, defined to be = qi'y^b as in 
Eq. (2.7), is described in HQET by 


= Cv^^v^qK + Cv^qi'j'iK -, 

i=l 


(2.26) 


^In NRQCD applications the relative velocity between the heavy quark and heavy anti-quark 
should not be confused with the velocity v introduced here. 



where is the HQET field, which satisfies Eq. (2.19) and whose dynamics are given 
by £hqet- The dimension-four operators are 


■VI = -v>^qp±hy, 

(2.27) 

■V2 = QillP^hy, 

(2.28) 

■V3 = 

(2.29) 

\v4 = +v^qp^hy, 

(2.30) 

■V5 = Qp±hlhy, 

(2.31) 

pQ = qiD^^hy. 

(2.32) 


Further dimension-four operators are again omitted, because they are linear combinations of 
those listed and others that vanish by the equations of motion. For example, q{iv ■ D)v^hv = 
— p)v^hy = Qy^ — rriqV^qhy, where the Dirac equation is used for the last step. 

The axial vector current has a completely analogous description, 

6 

A^ = CA^qiliTi^hy - CA^v^qi^shv - ^aiQm H-> (2-33) 

i=l 


where each operator is obtained from by replacing q with —gys. 

The short-distance coefficients of HQET depend on the heavy-quark mass m, as well as 
/i/m and mq/m, where /i is the the renormalization scale and fUq is the light quark mass. 
They are not explicitly needed in this paper, but it may be instructive to give the coefficients 
of the dimension-three terms through one-loop order, with niq = 0 (J = H, H) [27]: 

C'jii = 1 + (ih - 2 ) , (2.34) 

Cj^ = l + (7hln(mVh^) - 4 ) , (2.35) 

where the anomalous dimension 7 /^ = 3/2. The /i-independent part of Ca^^ and Ca^ given 
here assumes that the axial current is renormalized in a chirally symmetric way [38]. The 
coefficients of the dimension-four currents are 

= *<2. (2.36) 

s>3, (2.37) 


at the tree level, but 


all Bji become non-trivial when radiative corrections are included. 


C. HQET description of lattice gauge theory 

HQET provides a systematic way to separate the short distance 1/m from the scale Aqcd 
in heavy-light matrix elements, as long as the condition (2.16) holds. The formalism can also 
be applied to lattice gauge theory, again as long as condition (2.16) holds (and pa 1 ). 
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When lattice NRQCD is used for heavy-light systems, this is because £hqet is just the 
Symanzik LE£ for lattice NRQCD. When Wilson fermions are used for heavy quarks, one 
may also apply HQET, because they have the same particle content and heavy-quark sym¬ 
metries [4], In both cases bilinears of lattice fermions fields are introduced to approximate 
the continuum QCD currents. One held corresponds to the light quark, and the other to 
the heavy quark. An explicit construction, through order 1/m, is in Ref. [29] for lattice 
NRQCD, and a similar construction for Wilson fermions is in Sec. III. Lattice artifacts 
stemming from the light quark can be described as in Sec. II A, but lattice artifacts of the 
heavy quark should be lumped into the HQET short-distance coefficients. Some of the oper¬ 
ators needed to describe heavy-quark discretization effects do not appear in the usual HQET 
description of continuum QCD. Eor example, the dimension-seven operator h^Dfhv (writ¬ 
ten here in the rest frame) appears in to describe the breaking of rotational invariance 
on the lattice. Similarly, at and beyond dimension hve there are HQET current operators 
to describe violations of rotational symmetry in the lattice currents. Because of the high 
dimension, these effects lie beyond the scope of this paper, which concentrates on operators 
of leading and next-to-leading dimension. 

In this way, the preceding description of continuum QCD can be repeated for lattice 
gauge theory with the same logic and structure. Instead of Eq. (2.1), one introduces a 
relation like Eq. (2.18), 


-^lat — -CHgET, (2.38) 

where £iat is a lattice Lagrangian for NRQCD or Wilson quarks, and £hqet is an HQET 
Lagrangian with the same operators as in Eqs. (2.20) and (2.21), but modihed coefficients. 
In the dimension-four HQET Lagrangian one must now replace m with the heavy quark 
rest mass mi. The other coefficients will be denoted C*®'*. In particular, in the coefficient 
of the kinetic energy becomes 

^lat ^ 

^ 2m2 

If operator insertions of O 2 continue to be dehned by dimensional regularization with mini¬ 
mal subtraction, then both the rest mass mi and the kinetic mass m 2 generalize the pertur¬ 
bative pole mass. Like the usual pole mass, they are properties of the pole in the perturbative 
quark propagator [31], and they are infrared finite and gauge independent [39]. The lattice 
breaks Lorentz (or Euclidean) invariance, so reparametrization invariance no longer requires 
m 2 to be the same as mi. 

Similarly, a heavy-light lattice (axial) vector current (^ilt) can be described by 

r,i‘, = c'^yqK + c|?yn.y!„ - Y, + ■ ■ •, 

i=l 

At, = - E B'iQ% + ■ ■ ■, 

i=l 

but there are two important changes from Eq. (2.26). First, the light quarks (and gluons) 
are now also on the lattice, so they are described by their usual Symanzik LE£s. Second, the 


(2.40) 

(2.41) 
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short-distance coefficients of HQET are modified, because the lattice modihes the dynamics 
at short distances. The coefficients and now depend on the lattice spacing a, 

i.e., on ma, in addition to m, /r/m, and mq/m. A heavy-light lattice axial vector current 
has an analogous description. 

On the other hand, in Eqs. (2.26) and (2.40) the HQET operators are the same. As a 
rule, the ultraviolet regulator of an effective theory does not have to be the same as that of 
the underlying theory. (The standard Symanzik program works this way.) Thus, when de¬ 
scribing lattice gauge theory one is free to regulate HQET just as one would when describing 
continuum QCD. Moreover, since Eqs. (2.27)-(2.32) give a complete set of dimension-four 
HQET currents, the coefficients and contain short-distance effects from both the 

light and the heavy sectors. 

By comparing the HQET descriptions of lattice and continuum QCD, one can see how 
lattice matrix elements differ from their continuum counterparts. The continuum matrix 
element of n ■ V, for example, is 

{L\v ■ V\B) = -Cv,{L\qK\B^^^) - Bvi{L\v ■ - Bva{L\v ■ Qv4\B^^^) 

- I d*x{L\T02{x)qK\Bi^^y - CbCv^^ J d"a;(L|T 06(a;)gh,|Hf )* 

+ 0(AVmQ, (2.42) 

where L is any light hadronic state, including the vacuum. (The T*r-ed T product is dehned 
in Ref. [4]; this detail is unimportant here.) On the left-hand side B denotes a 6-flavored 
hadron, and on the right-hand side B^j^^ denotes the corresponding eigenstate of the leading 
effective Lagrangian Similarly, the lattice matrix element is [4] 

{L\v ■ V,,,\B) = - B^^\{L\v ■ Qyi|5®) - B^^l{L\v ■ Qy4|5®) 

J d‘^a:{L\T 02 {x)qK\Bi^'>r-C^i^C]^; J d^x{L\T OB{x)qK\Bi^'>r 

- j d^x{L\T qiaFq{x)qK\Bi^^y + 0{A^a^b{ma)). (2.43) 

Compared to Eq. (2.42), the short-distance coefficients are modified to depend on ma, there 
is an extra term from the Symanzik LE£ of the light quark, and the next power corrections 
can, in general, be multiplied by a (bounded) function of ma. The matrix elements on the 
right-hand sides are, however, identical, because in both cases they are dehned with 
describing the heavy quark and Dqcd describing the light quark (and gluons). 

Similar equations hold for matrix elements of Vj_ and VAiat, and for the axial vector 
current. If one multiplies the equations for the lattice matrix elements with 

Zj, = ( 2 - 44 ) 

Zj, = (2.45) 

and subtracts the result from the continuum equations, one hnds that the difference can be 
traced solely to the mismatch of the short-distance coefficients, or 
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<5C, = C]"* - C,, (2.46) 

5Bj, = - B,n, (2.47) 

where the normalization factors Zji are for i = 1, 4, and Zj^ for i = 2, 3, 5, 6. In 
Eqs. (2.46) and (2.47) a picture emerges, where heavy-quark lattice artifacts are isolated 
into SCi and SBji. Furthermore, the analysis presented here makes no explicit reference 
to any method for computing the short-distance coefficients, so it applies at every order in 
perturbation theory (in and, presumably, at a non-perturbative level as well. 

The matching factors Zj^^ and Zj^ play the following role, sketched in Fig. 1. In each 
case, the denominator converts a lattice-regulated scheme to a renormalized HQFT scheme, 
and the numerator converts the latter to a renormalized (continuum) QCD scheme. As long 
as the same HQFT scheme is used, HQFT drops out of the calculation of Zj^^ and Zj^. 
Moreover, changes in continuum renormalization conventions modify only the numerator, 
and changes in the lattice action or currents modify only the denominator. In a similar way, 
dependence on the HQFT renormalization scheme drops out when computing SCi and SBji. 

One can derive a connection between the matching coefficients of the HQFT and the 
Symanzik descriptions when ma -C 1 and m 3> p, so that both formalisms apply. With 
the Lagrangian, one applies HQFT to Fqs. (2.2)-(2.4) and identihes the short-distance 
coefficients with mi, and Cg*. Then one hnds, 

mib = mb + O(a^), 
m2b = mb + O(a^), 

= Zb- ^mbaK^.pC^.F, 

where the short-distance coefficient appears in the relation 

= -2C^.fOb. (2.51) 

At the tree level, C'f^F = 1- For the [axial] vector current, one inserts Fq. (2.26) [Fq. (2.33)] 
into Fq. (2.5) [Fq. (2.6)], neglects terms of order m'^af, and compares with Fq. (2.40) 
[Fq. (2.40)]. One also must match the tensor and pseudoscalar bilinears to HQFT at the 
dimension-three level, 

qia^’^b = - CT_q{v^iY± - (2.52) 

giysfe = (2.53) 


(2.48) 

(2.49) 

(2.50) 



FIG. 1. Diagram illustrating how the matching factors C, and Z 
gauge theory and QCD to HQET, and to each other. 


Cjc^a.t lattice 
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with short-distance coefficients Ct± and Cp. At the tree level, = Cp = 1. After 
carrying out these steps, one hnds that 

^V|| = Zy, (2-54) 

Zy^ = Zy^ (mg -h mb)aKyCT_/Cy^, (2.55) 

Zy^^By\ = Byi -\- aZyKyCT_-, (2.56) 

^V\\Byi = Byi + aZyRyCp^, f = 2,6, (2-57) 

^V||-By 3 = By^ — aZyKyCpp, (2.58) 

Zy^^By\ = Byi — aZyKyCT_, (2.59) 

Zy^^Bf, = By, - aZyKyiCpp - Ct_) (2.60) 

from matching the vector current, and 

= ^A, (2.61) 

^A|| = ^A^ + + mb)aKACp/CA^^, (2.62) 

ZA^B^2' = BM + 0{a^), i = l,2,5, (2.63) 

^A±B^2i = BAi + aZAKACp, f = 3, 6, (2.64) 

Za^B^aa — Baa ~ qZaKaCp (2.65) 


from matching the axial vector current. Of course, these relations hold only when describing 
the same lattice currents and and then only to order a^. Considering similar 
relations for the whole tower of higher-dimension operators, one sees 

hmC^* = Co, (2.66) 

a^O 

\imZj,B^^l = B,n. (2.67) 

a^O 

Eqs. (2.55)-(2.60) and (2.62)-(2.65) illustrate for the next-to-leading dimension operators 
how the limit is accelerated for standard 0{a) improvement, with K^.p, Ky, and Ka them¬ 
selves of order a. 

Equations (2.54)-(2.65) show that HQET matching connects smoothly to Symanzik 
matching in the limit where both apply. HQET matching is, therefore, a natural and 
attractive extension into the more practical region where ma is not very small. Continuum 
QCD still can be approximated well, but now order by order in the heavy-quark expansion. 

The remainder of this paper pursues this program in perturbation theory. One-loop 
corrections to the rest mass mi and the kinetic mass m 2 have been considered already in 
Ref. [31]. The one-loop correction to Cp would require a generalization of the calculation 
of K„.p [33] to incorporate the full mass dependence of the quark-gluon vertex. In this 
paper we focus on heavy-light currents. We construct lattice currents suitable for matching 
through order l/mg in the heavy quark expansion. We then calculate the matching factors 
Zj|| and Zj^ at the one-loop level, which are needed to £x the overall normalization of the 
heavy-light currents. Currents suitable for heavy-to-heavy transitions b ^ c are considered 
in a companion paper [6]. 
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III. LATTICE ACTION AND CURRENTS 


In this section our aim is to define heavy-light currents with Wilson fermions that are 
suited to the HQET matching formalism. Because Wilson fermions have the right particle 
content and obey the heavy-quark symmetries, the descriptive part of the formalism applies 
in any case. To use HQET to match lattice gauge theory to continuum QCD, however, we 
would like to ensure that SCi and SBji [cf. Eqs. (2.46) and (2.47)] remain bounded in the 
infinite-mass limit. Good behavior is attained by mimicking the structure of Eqs. (2.27)- 
(2.32), so that improvement terms are guaranteed to remain small. Then we would like to 
adjust free parameters in the currents so that 6Ci and SBji (approximately) vanish. We 
show how to do so in perturbation theory, obtaining B^f^ at the tree level and, in Sec. IV, 
the matching factors and Zj^ at the one-loop level. 

A suitable lattice Lagrangian was introduced in Ref. [5]. It is convenient to write the 
lattice Lagrangian £iat = Bq + The first term is 


To = -(mo + mocr)^(x)^(x) - |^7(a:) [(1 + 74)T>4 - (1 - 74)T>^iat] 

- C^(a^)7 ■ Di^ti’(x) + ^rXa^(x)Al^l'i/j(x). 


(3.1) 


The mass counterterm mocr is included here so that, by definition, mo = 0 for massless 
quarks. The covariant difference operators Di^t, and are defined in Ref. [5]. 

They carry the label “lat” to distinguish them from the continuum covariant derivatives in 
Secs. IIA and IIB. The symbol ip is reserved in this paper for lattice fermion fields. The 
temporal kinetic term is conventionally normalized, but the spatial kinetic term is multiplied 
with the coupling Q The coupling is, in the technical sense, redundant [5], but is included 
to solve the doubling problem [17]. 

For To the tree-level relations between its couplings and the coefficients in the Thqet 
are well known. By matching the kinetic energy, one finds (for v = 0) 

^lat[0] _ ^ _ _C_ I _/g 2^ 

^ 2 m^ 2 ^a moa{2 + moa) 2{l + moa)' 

At higher orders in perturbation theory, C 2 remains (for ?; = 0) the kinetic mass of the 
quark, which is expressed in terms of the self energy in Ref. [31]. 

To has cutoff artifacts, which are described by dimension-five and -higher operators in 
'^Sym (if mqtt 1) or Trqet (T 3> Aqcd)- The dimension-five effect can be reduced by 
adding 


Cb = facsCV’(a^)S ■ Bi^t{x)ip{x), (3.3) 

Ce = \acEC, ip{x)oL ■ E\^t{x)ip{x), (3.4) 


and suitably adjusting of cb and ce- The lattice chromomagnetic and chromoelectric fields, 
Slat and Slat, are those given in Ref. [5]. 

By matching the gluon-quark vertex, one finds 


^lat [0] 


1 

2mg^a 


moa(2-|-moo) 2(l-|-moa) 


(3.5) 
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Higher-order corrections to Cg* have not been obtained. By comparing Eqs. (3.2) and (3.5) 
one sees, however, that cb = rg + 0{g^) is needed to adjust Cg* to its continuum counter¬ 
part Cb = Zb I‘21712- 

The Euclidean action is S' = Special cases are the Wilson action [17], 

which sets = C = 1, cs = = 0; and the Sheikholeslami-Wohlert action [16], which 

sets Vg = ( = 1, cb = ce = csw- But to remove lattice artifacts for arbitrary masses, the 
couplings Tg, (, cb and ce must be taken to depend on moa [5]. Our analytical results for 
the integrands of Feynman diagrams, given in Appendix B, are for arbitrary choices of these 
couplings. Indeed, our expressions allow the heavy and light quarks to have different values 
of all couplings. 

Heavy-light currents are dehned in an essentially similar way. For convenience, first 
dehne a “rotated” field [18,5] 

Tg = [1 -h adi7 ■ Uiat] i^q, (3.6) 

where ijjq is the held in Co of havor q, and Z^iat is again the symmetric covariant difference 
operator. Simple bilinears with the right quantum numbers are 

(3.7) 

AO = Tgiy^ysTfc. (3.8) 


The subscript “0” implies that, as with Cq., some improvement is desired. To ensure a 
good large-ma limit, one should pattern the improved current after the right-hand side of 


Eq. (2.40). Thus, we take 

( 3 . 9 ) 

i=l 

K, ^ At; - ( 3 , 10 ) 

1 = 1 

where the bj; are adjustable, and the dimension-four lattice operators are 

Qvi = (3.11) 

Qv2 = ( 3 - 12 ) 

Qv3 = (3.13) 

Qvi = (3.14) 

Qvb = ^qPuJll'^b, (3.15) 

Qv& = (3-16) 


and each lattice operator is obtained from Qy^ by replacing ijig with — V’g 75 . Lattice 
quark helds do not satisfy Eq. (2.17), so ^ appears explicitly. In practice, one uses the rest 
frame here, v = (f,0), as in Eq. (3.6). An analogous construction for lattice NRQCD has 
been given by Morningstar and Shigemitsu [29]. 

It is worthwhile to emphasize the difference between Eqs. (2.40) and (3.9). Equa¬ 
tion (2.40) is a general HQET description of any heavy-light lattice current. Equation (3.9) is 
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a definition of a specific lattice current, namely the one used in this paper (and in calculations 
of fs and other hadronic matrix elements). In the same vein, the Qji in Eqs. (2.27)-(2.32) 
are HQET operators, whereas the Qji in Eqs. (3.11)-(3.16) are lattice operators. Finally, 
the coefficients are the output of a matching calculation: they depend on the bji, which 
must be adjusted to make SBji vanish. 

To illustrate, let us consider the calculation of the coefficients B]jl at the tree level. One 
computes on-shell matrix elements such as {q\Jia.t\b) and {0\Jiat\qb) in lattice gauge theory 
and compares them to the corresponding matrix elements in HQET. Then one hnds 


^lat[0] 

J\\ 


^ Jl_ ’ 


B 


latIO] 

Ji 


B 


latIO] 

Ji 


e 

e 


(m™+m™)a/2 f _J: _, ^[0]\ 


i < 2 


(3.17) 

(3.18) 

(3.19) 


where 


= ln(l -|- moa) 

and, for our lattice Lagrangian and currents, 

1 _ ((1+moa) 

2m[j°^a moa(2 + moa) 


(3.20) 


(3.21) 


Since (continuum QCD’s) Cj = 1 there already is a non-trivial matching factor at the tree 
level relating the lattice and continuum currents, Zj^j = Zjj^ = e(™'i 9 +™i 6 )“/ 2 _ 

After comparing Eqs. (3.18)-(3.19) with Eqs. (2.36)-(2.37), one sees that one can take 
6jj = 0 for all six operators, if di is adjusted correctly. At the tree level, the way to adjust 
di is to set equal to the (tree-level) heavy-quark mass. In the effective Lagrangian 
there are two quark masses, the rest mass mi and the kinetic mass m 2 - The former has no 
effect on matrix elements (and a trivial, additive effect on the mass spectrum). As discussed 
above, heavy-quark cutoff effects in matrix elements are reduced if = C 2 , which means 
one should identify the continuum quark mass with the kinetic mass. Thus, one should set 
m^o^ = m ^2 \ which is obtained if one adjusts 

d, = , (3.22) 

moa{2 + moa) 2{1 + moa) 

The same rotation also improves heavy-heavy currents at the tree level. 

Beyond the tree level, it is convenient to define di for the spatial component of the 
degenerate-mass, heavy-heavy current [6]. Then the corrections heavy-heavy current analo¬ 
gous to Qv 2 and Qvb would be superfluous, but for unequal masses they are still required. 

For equal mass currents it is possible to compute Zy^^ nonperturbatively for all masses mi,. 
One may therefore prefer to write [13-15] 



, / Zyuu Z\rbb P Jut 

V bi bi 


ub 


(3.23) 
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and compute only the factor pjub in perturbation theory. To calculate the pre-factor Zybb 

appearing in Eq. (3.23), one must have a massive quark in the hnal state. The dehnition of 
the heavy-heavy matching factor is given in our companion paper [6], along with a calculation 
of its one-loop level contribution. We give the results for heavy-light ^ in Sec. IV. 

For a light quark, with m^a 1, the right hand side of Eq. (3.22) vanishes linearly 
in moa. Therefore, 1 — mzq/m 2 q is 0 ( 171 ^ 0 ^), and the distinction between m^q and m 2 q is 
negligible. For this reason, and to simplify calculation, we set niq = 0. Then Eq. (3.22) 
implies di = 0 for the light quark. 


IV. ONE-LOOP RESULTS 

In this section we present results for the matching factors at the one-loop level in pertur¬ 
bation theory. The one-loop contributions are known for the Wilson [11] and Sheikholeslami- 
Wohlert (SW) actions [12]. Both these works omit the rotation term in the current [18,5], 
which is needed to obtain l/ms correctly. In this section we complete the work started in 
Ref. [12] and report results with the clover term and with the rotation. For comparison we 
also present our results without the rotation, both with and without the clover term. 

The computer code for generating these results is freely available [30]. 

The matching factors Zj {J = V\\, V±, A\\, and A±) are simply the ratios of the lattice 
and continuum radiative corrections: 


71/24 71/2 

cont 

7I/2A 7I/2 

^ 2 h 

lat 


(4.1) 


where Z 2 h and Z 21 are wave-function renormalization factors of the heavy and light quarks, 
and the vertex function Aj is the sum of one-particle irreducible three-point diagrams, in 
which one point comes from the current J and the other two from the external quark states. 

The expression relating Z 2 to the lattice self energy, for all masses and gauge couplings, 
can be found in Ref. [31]. Its dominant mass dependence is 


Z 2 (X e 


—mia 


( 4 . 2 ) 


where rrti is the all-orders rest mass (of the heavy quark). This mass dependence is not 
present in the vertex function or the continuum part of Eq. (4.1). Consequently, we write 

rni 

e-™. (4,3) 

1=1 

SO that the Zjj, are only mildly mass dependent. (A slightly different convention was used in 
presenting results for Z 2 in Ref. [31].) By construction, this mass dependence in pj^ cancels 
out in a gauge-invariant, all orders way. So, we write 

CXD 

PJr = 1 + X^PoVl/r- (4-4) 

l=l 
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This rest of this section is split into two subsections. In the first, we present our results 
for the full mass dependence of and f)j\ In the second, we discuss the related calculation 
of the Brodsky-Lepage-Mackenzie scale q*. In both cases, we discuss fully a range of checks 
on our calculations. 


A. and 

The combinations of wave-function and vertex renormalization in Zj are gauge invari¬ 
ant and ultraviolet and infrared finite. For vanishing light quark mass there is a collinear 
divergence (which can be regulated by an infinitesimally small mass), but it is common to 
lattice and continuum functions. In the desired ratio (4.1), the divergence cancels, and the 
result is independent of the scheme for regulating the collinear singularity. For large ma a 
remnant of this cancellation appears. The lattice theory approaches its static limit, where 
its ultraviolet behavior is non-logarithmic. But the region of momentum a~^ < q < m in 
the continuum diagrams generates logarithms. At the one-loop level one must find 3 In(ma), 
with the same anomalous dimension as in Eqs. (2.34) and (2.35). At higher loops the usual 
polynomial in In(ma) will arise. 

We have calculated the one-loop Feynman diagrams for the action specified in Eqs. (3.1)- 
(3.4), with arbitrary mo, Ts, (C, cb, and Ce for the incoming heavy quark, mg = 0, r(, 
C', c^, and for the outgoing light quark. The needed Feynman rules are in Ref. [31], 
apart from three new rules for the current itself, which are in Appendix A. As shown in 
Appendix B, we have found a simple way to incorporate the rotation into the Dirac algebra. 
The resulting analytical expressions are surprisingly compact, and they are given explicitly 
in Appendices B and C. 

We have evaluated these expressions for = C = 1 and ce = cb = csw- Thus, the 
numerical results correspond to the SW action (csw = 1) and to the Wilson action (csw = 0). 
Figure 2 plots the full mass dependence of the matching factors for the vector current, 
(a) (b) Zy^ , (c) pyi, and (d) pvy • These numerical results are for the SW action with 

rotation (solid lines) and also for the SW and Wilson actions without the rotation (dotted 
lines). Figure 3 plots the full mass dependence of the matching factors for the axial vector 
current, (a) Z^||, (b) Za^, (c) pAy, and (d) p^^. These and the following figures are plotted 

against mf^a because this variable conveniently covers the whole mass range: for small mass 
mi ~ m 2 , and for large mass rriia ~ lnm 2 a. 

We have carried out several checks on our calculations. In each case, identical numerical 
results have been obtained with two or more completely independent programs. The results 
for agree with those previously obtained, for csw = 0 [11] and for csw = 1, di = 0 [12]. 
We have also reproduced limiting cases, as we briefly discuss below. 

For ma = 0 our calculation reduces to the usual matching calculation for massless quarks. 


We find (with Ce = 4/3) 




^V|l - - j 

[ -0.129423(6), 

[ -0.174073(7), 

Csw — 1; 

Csw = 0, 

(4.5) 

z'a, = Za[ = 1 

[ -0.116450(5), 

[ -0.133365(5), 

Csw = 1; 

Csw = 0, 

(4.6) 
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FIG. 2. Full mass dependence of the one-loop coefficients of the matching factors of the vector 


current (a) Zy^, (b) Zy^, (c) and (d) py|. Filled (open) symbols denote the SW (Wilson) 
action; solid (dotted) lines connecting squares (circles) indicate the rotation is included (omitted). 


in excellent agreement with previous work for csw = 1 [40-42] and csw = 0 [42-44]. (Refer¬ 
ence [42] gives precise results as a polynomial in cswO 

As the mass tends to infinity, these actions and currents all lead, up to an unphysical 
factor, to the same vertices and quark propagator—a Wilson line. Perturbative corrections 
to the vertex functions must respect this universal static limit, and, therefore, they must 
tend to a universal value. As ma oo, one expects the Z factors for a massive quark to 
approach those for the static limit, namely 

Zf = \lh \n{m2af + , (4.7) 

where the constant Zj^ depends on the current J and on csw (of fhe light quark). Since 
ln(m 2 a) ~ rriia in this region one expects the linear behavior seen in Figs. 2 and 3. The 
static limit is also shown in Figs. 2 and 3 with 


[1] ^ I —10.248, Csw = 1, 

V|| - j _7.929, Csw = 0, 

^[1] ^ f —12.248, Csw = 1, 
I —9.929, Csw = 0, 


(4.8) 

(4.9) 


19 








FIG. 3. Full mass dependence of the one-loop coefficients of the matching factors of the axial 


vector current (a) (b) zj^[, (c) , and (d) 

[1] ^ f —18.414, Csw = 1, 

\ -24.379, Csw = 0, 

[1] ^ f —16.414, Csw = 1, 

^11 “ \ -22.379, Csw = 0. 

We have obtained these constants ourselves [55]. They agree with previous (less precise) 
results for csw = 1 [45] and csw = 0 [27]. As one can see from looking at Figs. 2 and 3, the 
static result is a good approximation for mf'a > 5 or, equivalently, moa ~ 77120 > 150. 

Some of the points at the highest masses have large error and he nearly one a off the 

curve. The origin of this behavior is that the lattice and continuum integrals are dominated 
by different momenta: the continuum integral is dominated by the region k ~ m 2 ^ a~^, 
whereas the lattice integral is dominated by the region k ~ a~^. This mass region is not of 
much practical interest, since here one has an essentially static quark. 

Equations (2.54)-(2.65) allow us to check the small (heavy-quark) mass limit against 
the work of Sint and Weisz [46]. In our conventions the matching factors Zy and Za are 
functions of gauge coupling and quark mass. Thus, 

Zy{mqa,mi,a) = Zy 1 + \{mq +mi,)aby , (4.12) 

ZAimqtt, mta) = Za I + \ {mq + mh)abA , (4-13) 


(4.10) 

(4.11) 
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where, on the right-hand side, we adopt the notation of Refs. [9,10,46], and the Zs and bs 
do not depend on mass. Here only the mass dependence is displayed; all quantities depend 
also on the gauge coupling. 

If we omit the rotation, our currents and those considered by Sint and Weisz coincide, 
apart from one-loop counterterms. Thus, in one-loop calculations the slopes of our mass- 
dependent matching factors must agree with them. Setting rriq = 0, and using Eqs. (2.54)- 
(2.65), 


dZ^ 


[ 1 ] 

Vii 


16 


dm 


dm 

dZ 


dm 


dZ 


dm 


— 1 


16 

[ 1 ] 


16 

[ 1 ] 

Ax 


ifth] _ ^[1] 
2'^V 1 


ifth] _ ^[ 1 ] 
2^A 1 
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2^A ■ 


( 4 . 14 ) 

( 4 . 15 ) 

( 4 . 16 ) 

( 4 . 17 ) 


To extract these slopes, we form a combination of integrands with three different (small) 
values of mbO, yielding and R'j' up to O(mfea)^. In this way we hnd (for csw = 1) 


feW = X 0.114929(10) = 0.153239(14) (4.18) 

vs. Cf X 0.11492(4) [46], 

5 W = CfX 0.114142(10) = 0.152189(14) (4.19) 

vs. Cf X 0.11414(4) [46], 

= Cf X 0.0122499(6) = 0.016332(7) (4.20) 

vs. Cf X 0.01225(1) [46], 

= CfX 0.0056806(11) = 0.0075741(15) (4.21) 

vs. Cf X 0.005680(2) [46], 


which agrees perfectly with Ref. [46]. These results have also been checked by Taniguchi 
and Ukawa [47]. We also obtain 

5 W - 5 W = CfX 0.0007833(11) = 0.0010444(16) (4.22) 


by subtracting the integrands hrst, and then integrating. In taking the difference, large 
contributions from the self energy cancel, but, even so, the near equality of b^y and b^^ is a 
bit astonishing. Comparing the slopes of Figs. 2(a) and 3(b) one sees that — b^^ for the 
Wilson action is not so small. 

Although these checks are reassuring, the main result of this section is to obtain the full 
mass dependence of the matching factors. The results at intermediate mass, with mia < 3 
or, equivalently, moa < 1.5, are needed for realistic calculations of B meson properties. This 
region is neither particularly close to the massless limit, nor to the logarithmic behavior of 
the static limit. 
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B. BLM scales q 


It is well-known that perturbation theory in the bare coupling g^ll/a) converges poorly. 
Therefore, we calculate the ingredients needed to determine the Brodsky-Lepage-Mackenzie 
(BLM) scale [48,49]. For a coupling in scheme S, we denote the BLM expansion parame¬ 
ter gsi^s)- The BLM scale is given by 


In^ggaY 


( 1 ) J d^k \n{kaf f{ka) 
^ f d'^kf{ka) 


(4.23) 


where k is the gluon momentum, and f{k) is the integrand of the quantity of interest, e.g., 
/ d'^k f{k) = Zj'K The constant is the /5o-dependent part of the one-loop conversion 
from the arbitrary scheme S to the ‘V scheme”, namely 


9s{(l) 


9U<1) 


+ P„bf+bf+0(g^), 


(4.24) 


where for Uf light quarks /do = 11 — 2n//3, and 6® is independent of n/. The V^-scheme 
coupling gy{q) is dehned so that the Fourier transform of the heavy-quark potential reads 
V{q) = —CF9vi.9)/<f- Equation (4.23) shows that the dehnitions of q* in Refs. [48] and [49] 
are identical in the V scheme. 

For our matching factors it is straightforward to weight the integrands with ln(A;a)^ to 
obtain 

* T’ll] 

In(gya)^ = ^-fi] ’ (4-25) 


because the integration over d^k has no divergences. The denominators are the one-loop 
coefficients given above, and the numerators are presented now. 

Figure 4 plots the full mass dependence of the numerators for the vector current, (a) , 

(b)* Zy^, (c) and (d) As before, these numerical results are for the SW action with 
rotation (solid lines) and also for the SW and Wilson actions without the rotation (dotted 
lines). Figure 5 plots the full mass dependence of the numerator of Eq. (4.25) for the axial 
vector current, (a) (b) (c) *p! 4 ||, and (d) *We have carried out several 

checks on our calculations. Once again, identical numerical results have been obtained with 
two or more completely independent programs. Also, at mi,a = 0 we reproduce the results, 
for the Wilson action, of Ref. [44]. 

For *z'f^ and the limit of large ma also has distinctive features. In that case 


* 


zf 


Cf 


l'yh\n‘^im2af+ 'y'f^\n{m2ay+ *zj , 


(4.26) 


where is related to the two-loop anomalous dimension. A similar expression holds for 
*p^j\ with a different constant. Note that—in both cases—the one-loop anomalous dimension 
appears multiplying ln^(m 2 a). The growth expected from Eq. (4.26) is seen in Figs. 4 and 5. 
As a consequence, one hnds q*a oc as ma oo. Square root behavior is typical of 

cases with an anomalous dimension. 
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FIG. 4. Full mass dependence of the estimated terms of the matching factors of the vector 
current (a) (b) (c) and (d) . 


For the Z factors, the resulting values for q*a are relatively constant in the “low mass” 
region, q*a ~ 2.7-2.9. Figure 6 shows how q*a depends on the heavy quark mass in the 
region mia < 2, which is the one most relevant to calculations of decay constants and 
form factors. At larger masses goes through zero, at which point the original BLM 

prescription breaks down. A prescription for q* in this case is given in Ref. [50]. For the 
Wilson action the zero in Zy^ is at a smaller than usual mass [see Fig. 2(a)], which explains 

its behavior for the BLM q*a seen in Fig. 6(a). For the p factors the denominator is 
small over most of the interesting region, as seen in Figs. 2(c)-(d) and 3(c)-(d). 

It is also interesting to see how q* changes under tadpole improvement. If one introduces 
the tadpole-improved matching factors 

Zj = Zj/uo, (4.27) 

where the mean link uq is any tadpole-dominated short-distance quantity, the arguments of 
Ref. [49] suggest that the perturbative series for Zj has smaller coefficients. In analogy with 
Eq. (4.1) we write 

rni 

^/^Zj = l + Y,gl^z'f, (4.28) 

i=i 


where 
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FIG. 5. Full mass dependence of the estimated terms of the matching factors of the axial 
vector current (a) (b) *Z^[, (c) and (d) 


rrii^a = ln[l + rrioa/uo] (4.29) 

is the tadpole-improved rest mass. Then 

4 “' = 4 ' - ) (i + — ) “™. (4'30) 

2 V 1 -b m^aj 

and because z'j^ < 0 and Mg' < 0 one sees that the one-loop coefficients are reduced. 
Similarly, for computing the BLM scale 

To illustrate, we take mq from the average plaquette, so Mg' = —Cf/ 16 and *Mg^ = 
—0.204049(1). Figure 7 shows that, as a rule, q* is significantly reduced, which means 
that tadpole improvement has removed some of the most ultraviolet contributions. With a 
lower scale, the coupling gy{q*) becomes a bit larger with tadpole improvement. For 

and Za^i however, the denominator z'j^ already vanishes for mia ~ 1.5-2.0, leading to 
rapid growth in the BLM q* for the Wilson action, and a zero in the BLM q* for the SW 
action. One should again dehne q* in a more robust way [50]. Another choice for the mean 
held is Mg = S/tcrit- It gives coefficients and BLM scales that lie between the unimproved 
and tadpole-improved cases [51]. 
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FIG. 6. Full mass dependence of the BLM scale g*, for (a) Zv||, (b) (c) Za^^^-, and (d) Za^- 


Our method also allows us to obtain the BLM scale for the improvement coefficients in 
the Alpha program. Then we are in a position to compare BLM perturbation theory with 
non-perturbative determinations of these coefficients. We will give these results for q* and 
the mentioned comparison in another publication [52]. 


V. CONCLUSIONS 

In this paper we have set up a matching procedure, based on HQET, for heavy-light 
currents. It is valid for all ma, where m is the heavy quark’s mass and a is the lattice spacing, 
and to all orders in the gauge coupling. It could be applied to lattice NRQCD, although here 
it is applied to Wilson fermions. In the latter case, HQET matching agrees with Symanzik 
matching when ma 1. In this way, HQET matching is a natural and attractive extension 
into the regime ma ^ 1, which is needed for heavy-quark phenomenology. 

Our one-loop results for the SW action are of immediate value for lattice calculations 
of fs and of form factors for the semi-leptonic decay B ttIu. Indeed, our earlier one- 
loop results [12] (which omitted the “rotation” terms in the current) were used for fs 
in Refs. [19-23], and our results were used for semi-leptonic form factors in Ref. [15]. In 
particular, we have obtained the BLM scale q* for the matching factors, which should reduce 
the uncertainty of one-loop calculations. Similarly, computing part of the normalization 
factor, namely ^ZymZybb, non-perturbatively reduces the normalization uncertainty even 
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FIG. 7. Tadpole-improved q* for (a) Zvj|, (b) (c) Za^^-, and (d) Za^-, with Fig. 6 in gray. 


further [13-15]. (The heavy-heavy normalization factor Zyhh is dehned in our companion 
paper for heavy-heavy currents [6].) 

An outstanding problem at this time is the one-loop calculation of the coefficients 
of the dimension-four terms in the HQET description. A calculation of these coefficients, 
and the subsequent adjustment of the parameters bji in the lattice currents, would elimi¬ 
nate uncertainties of order agA/m and a^Aa in (future) calculations of heavy-quark matrix 
elements. The algebra quickly becomes voluminous, making this problem well-suited to 
automated techniques [53]. 
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APPENDIX A: FEYNMAN RULES 


The needed propagators and vertices for quark-gluon interactions are given in Ref. [31]. 
Here we give the additional Feynman rules induced by the rotation term of the heavy quark. 
The additional rules are easy to derive by expressing the covariant difference operator as [54] 

DL = [T*, - r-„] /(2o), (Al) 

where 

T±, = (A2) 

and t-i -^/2 translates fields to its right by one-half lattice spacing in the ±p direction. 

There are three rules to give, with 0, 1, and 2 gluons emerging from the vertex. Let the 
Dirac matrix of the current be T. Then, 


0-gluon = F 1 + sin(pr) , 

V 

(A3) 

1-gluon = gof^ di Fy* cos(p -h ^k)i, 

(A4) 

2-gluon = igl\{t‘",t’’}dij di F 7 jsin(p-h \k + |£)i. 

(A5) 


where momentum p is quark momentum flowing into vertex, and k and i are gluon mo¬ 
mentum flowing into vertex. As in Ref. [31], the matrices are anti-Hermitian, i.e., 


= exp (^go^A 




and trt“t^ = 


.l^ah 
2 " • 


APPENDIX B: DIRAC ALGEBRA 


To compute the vertex function, there are four diagrams to consider, depicted in Fig. 8: 
the usual vertex diagram (with the rotation inside). Fig. 8(a); two diagrams with the gluon 
connected to the incoming rotation. Fig. 8(b) and (c); and a tadpole diagram connected to 
the incoming rotation [using rule (A5)], Fig. 8(d). The tadpole diagram. Fig. 8(d), vanishes 
for zero external three-momentum, because i = —k and p* = 0. 

For each non-vanishing diagram. Figs. 8(a-c), dehne the integral 


T(a,b,c) 

-'r 



d k 1 (a,b,c) 

(^F ^ ’ 


(Bl) 




(b) 


(d) 


FIG. 8. Feynman diagrams for calculating the vertex function. The • on each side of the (8> 
indicates the rotation. 
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J 

r 

ST 

Vi 

74 

-1 

A4 

7475 

+1 


7i 

1 

3 


7i75 

+5 

TABLE I. 

The factor sp; defined by \ Yl,rlr^lr 

= stT. 


where k is the momentum of the gluon in the loop, and = 2sin(|A;^). Let the incoming 
massive quark have couplings mo, r-g, C) cs, and ce, and external momentum p. Similarly, 
let the outgoing massless quark have couplings mg = 0, r^, C', 4, and Cg, and external 
momentum p'. The internal quark lines carry momentum p + k in and p' + k out. The 
integrals I are obtained directly from the loop diagrams. Then 



( 7(1] _ 7LJ-J I yV _ 

\^'^2hcont ^2h\a.t ' ^21 cont 


dl] 


dl] 


- z. 


[ 1 ] 


2/ lat 


t+1: h 

d 


d 

r cont 



(B2) 


from Eq. (4.1). The relation between the current J and its Dirac matrix T is contained in 
Table I. The expression relating to lattice self-energy functions is in Ref. [31]. 

The most onerous task in evaluating the diagrams is the manipulation of the Dirac 
matrices. A convenient method is to treat each quark line separately, starting from the 
initial- or hnal-state spinor. Then the spinor, the propagator, and the vertices can be 
written out in 2 x 2 block diagonal form, with Pauli matrices appearing in the blocks. Once 
the Feynman rules are as complicated as in the present calculation, it is easier to manipulate 
2x2 matrices of Pauli matrices than to manipulate Dirac matrices. A special advantage of 
this organization is that the rotation bracket in Eq. (A3) merely “rotates” the rest of the 
leg. We also obtain in this way, with much less effort than in Ref. [31]. 

A further advantage is that the vertex corrections can be expressed compactly. The 
diagram with a gluon going from the incoming leg to the rotation. Fig. 8 (b), is 


n-ik) 




(3 





(B3) 


where Sr = sin/c^, and the functions D, L, and Kr are given in Appendix C. The diagram 
with a gluon going from the outgoing leg to the rotation. Fig. 8 (c), is 


T'ic) 


srcii-d 


+ \('Y.KS, 


(B4) 


where the functions D', L', and iP' are given in Appendix C, and sr is given in Table I. 
The unbarred function L (barred function L) is for P = 74 and 7^-75 (P = 7 ^ and 7475 ). 

The vertex diagram. Fig. 8 (a), is complicated. We hnd /DD\ with numerator 


nP = (±) (u\n[Uo] - srL'on[L^]S^) - CC'Xr. (B5) 


The upper sign and unbarred functions (lower sign and barred functions) are for P = 74 and 
7 j 75 (P = 7 j and 7475 ). The part Xr comes from spatial gluon exchange: 
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Xr = -sr(3 - ifc') + 4(3 - \k^)V'n[V]S^ 

+i (v'n[u] - srL'n[c]) Y.k.sI 

r 

+i (u'n[V] - »rC'7J|il) E KSl (B6) 

r 

+i (u'n[U] - srS^C'TJICl) E KkM 

r 

where the last term is absent for V 4 and A 4 (he., when Sp = 1). The rotation enters in the 
“rotated” functions 


7^[Uo] = Uo + dl5%, 

(B7) 

'T^[Lo] = Lq — diU^, 

(B8) 

n[u] = U + diS^C, 

(B9) 

7^[C] = C - d4U, 

(BIO) 

n[V] = U + diL, 

(Bll) 

7^[L] =L- diS^V. 

(B12) 


Although the vertex diagram is not easy to write down, the rotation modihes it in a fairly 
simple way, when using the 2x2 Pauli matrix method described above. 

We have verihed that these expressions are correct by completely independent calculation 
with more common methods for the Dirac algebra. 


APPENDIX C: USEFUL FUNCTIONS 

In this appendix we list the functions appearing in Appendix B for the action and currents 
given in Sec. III. First, let 


/i = 1 + mo + (Cl) 

ju' = 1 + m'o + . (C2) 

from now on a prime means to replace incoming couplings and momenta with corresponding 
outgoing couplings and momenta. 

When the quark propagator is rationalized it has the denominator 

D = 1 — 2/j.cos(k4 + + /j.^ + (C3) 

where = ln(l + mo). 

In this calculation, the heavy quark has zero three-momentum, so its spinor consists only 
of upper components. Depending on the matrix F the heavy quark couples either to the 
upper or lower components of the light quark. With the upper components the unbarred 
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functions arise, and with the lower components (of the light quark) the barred functions 
arise. 

To express the useful functions compactly, it is convenient to introduce hrst 


^ (C4) 

t/= ^ - (C5) 

because these combinations appear in the other functions. Then 

Uo = _ 1^2^^ cos{lh)S^, (C6) 

Lo = C + \ce cos(|A;4)f7j , (C7) 

V = Q 1 + |c£;sin(/c4) +\cbU, (C8) 

L =—U 1 + |c£; sin(A;4) +^cbCS‘^: (C9) 

Kr = rs- Cb cos^i^kr) = {Ts - Cb) + \cBkl, (CIO) 

and 

7/o = i/e—- 1(2^^ cos(iA;4)5^ (Cll) 

Lo = C ^ cos(iA;4)17] , (C12) 

C = C 1 — fc£;sin(A;4) + ^cbU (C13) 

L =—U 1 — Ics sin(A;4) +\cbC,S‘^. (C14) 

The barred functions are obtained from unbarred counterparts by putting k —k and 


so there is no need to introduce Kr = Kr. In the present calculation the 
barred functions arise only for the outgoing massless quark, for which ttIq = = 0. 
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